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Abstract 

In this work two dimensional simulations of ellipses models and realistic 
models of trabecular bone microstructure tissue are performed. The results 
presented here intend to show that microstructure is an important factor in 
the effective dielectric properties of trabecular bone. The main idea is to 
analyze the possibility of using the dielectric behaviour of bone tissue as a 
measurement of bone health. The frequency used was 950MHz. It was found 
that dielectric properties can be used as a estimate of degree of anisotropy, 
particularly, conductivity shown better results. For representing bone tissue 
effective models were also tested. Ellipses models with aspect ratio a/b = 
1.5 shown relative good agreement. As a final conclusion according to the re- 
sults presented here, anisotropy must be taken into account when measuring 
trabecular bone tissue with microwave imaging. 
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1. Introduction 

Osteoporosis is usually diagnosed by one parameter: bone mineral density 
(BMD). Unfortunately, it does not give direct tissue biomechanical informa- 
tion. There is a continuos search of several alternative methods for in vivo 
estimation of bone mechanical properties. In this regard, ultrasound seems 
to be the most advanced toolpQ. Other approaches, like that developed by 
Kanis and co-workers [2], intend to integrate the risks associated with clinical 
risk factors and BMD at the femoral neck. 

Measurement of dielectric properties of bones has been recently received some 
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attention due to its potential application as an alternative diagnostic method 
for evaluating bone health (3j HI El El [7J [E] . Low frequency in vitro measure- 
ments on a particular sample (human cadavers, age = 58 ± 20 years) have 
shown that linear combination of two relative permittivities in frequency 
ranges (50Hz-lkHz and 100kHz-5MHz), is able to predict 73% of the tra- 
becular bone volume fraction (ratio of bone volume to total tissue volume 
BV/TV) [9]. Clearly, at the state of the art, low frequency methods are 
only applicable to special cases of open surgery. Recently, the first relative 
permittivity and conductivity in vivo bone images have been obtained by 
microwave tomography [TUl [TT] . Furthermore, in vitro microwave studies 
suggest a negative correlation between bone mineral density and permittiv- 
ity / conductivity [T2l H5] . This agrees with the results of reference |14j . 
There, it is hypothesized that mineralized bones have lower relative permit- 
tivity mainly due to a lower water content. 

It is well known that not only mineral content but also microstructure defines 
the mechanical and dielectric properties of bone tissue [9] . There exist sev- 
eral parameters to characterize microstructure of bone, among others: BV 
/ TV, mean intercept length (MIL), trabecular thickness |1 5] . In reference 
[3] a theoretical approach for porosity estimation is developed. The authors 
use effective dielectric properties to obtain the spectral function of the bone, 
which is close related with its porosity. 

Age and pathological conditions produce selective bone tissue resorbtion [TB] 
(pag 326), which results in thinner horizontal trabeculae. Efforts have done 
in order to detect this anisotropy by using ultrasound techniques [T71 [18] . In 
those papers, the author uses fast and slow waves, which can propagate in 
the direction of the major trabecular orientation. Regarding this, trabecular 
bone was found to be electrically anisotropic. The relationships between the 
dielectric properties at low frequency (< 10MHz) in different directions have 
been early studied by Saha et. al [19]. These variations could be interpreted 
by trabecular orientation, but there is no strong evidence about it yet. How- 
ever, at low frequency (< 10MHz), no linear correlation between degree of 
anisotropy and dielectric or electrical parameters was found by Sierpowska 
et al. [9]. 

In electromagnetic simulation and theory of hierarchical materials, homoge- 
nization techniques are extensively used [20]. They give information of mi- 
crostructure [21] and anisotropy [22] . This paper intend to give a step toward 
a microstructure modeling of trabecular bone tissue. We study effective di- 
electric approaches and their application to realistic electromagnetic models 
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Figure 1: Simulation box. 

based on scaning electron microscopy images of bone tissue. The main idea 
is to study the anisotropy of trabecular microstructure and its relationships 
with the dielectric effective parameters using simulation techniques. The two 
dimensional (2D) models used here are numericaly solved by finite difference 
time domain (FDTD) [231 L21] • We will try to answer one important question: 
If we are interested in measuring trabecular bone at microwave range, is it 
confident to consider isotropic effective dielectric properties? At the same 
time, the answer can contribute to the development of a tool for microstruc- 
ture characterization of bone tissue by effective models which take account 
anisotropy. 



2. Methods 

Effective relative permittivity (e e s) calculus in 2D is study by observing 
the reflection of a slice of relative dielectric constant e e with inclusions of 
relative permittivity Figure [T] shows the simulation box. Two shapes of 
inclusion are tested: circular and elliptical. For validation purposes, different 
values of e e (<r e ) and E\ (a;) are tested (where a x is the conductivity of medium 
x). In all the cases overlaping is allowed. For simulation of realistic models of 
trabecular bone, we will also consider it as two phase material: the trabeculae 
as inclusions and bone marrow as background. Values of permittivity and 
conductivity are obtained from [25J, which is based on [26] . 

2.1. Numerical models 

Maxwell equations are numericaly solved using finite difference time do- 
main (FDTD) [21]. The implementation of FDTD algorithm is performed 
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using meep[2"5]. The slice is a square with variable thickness (d). The entire 
simulation box has variable length and it is adjusted large enough in order 
to reach three times the thickness of the sample at each side. A temporal 
gaussian pulse polarized in y direction (E y ) with central frequency equal to 
950MHz is used as source, located far enough from the interface air-slice. 
The resolution of the simulation box is 350/mi (or smaller) and the Courant 
number is S — 0.5. In x direction boundary condition is a perfectly matched 
load (PML). 

In order to obtain the e e ff the reflection coefficient of the medium with in- 
clusions or bone (T^u;)) is compared to an homogeneous block T(u,e e s)) 
in the frequency range of interest. Minimization of the absolute difference 
(|ri(o;) — r(u,e e fi)\) is performed using a Newton conjugate-gradient algo- 
rithm (implemented in Scipy 0.9.0 [27J). 

Scripts ad hoc in Python 2.7.2 |28j were developed to process the simulation 
data. In order to study different number of inclusions and also randomly 
locate circles and ellipses, ad hoc python routines were used as well. 

2.2. Effective permittivity models 

There are several models which represent effective dielectric properties 



of two-phase materials, Appendix A shows some of them. In these models, 



effective values depend on the fraction occuped by the inclusions {4>\) and 
on the permittivities (conductivities) e e (cr e ) and S\ (a{). They also depend 
on shape and anisotropy. Those models which take into account shape and 
anisotropy are of particular interest for this work [22], specifically elliptical 
inclusions will be studied here. The general form used for Mejdoubi and 
Brusseau [22] is: 

£ C R = e e -f(^, 0i, A^j « e e (1 + «0i + /30f) (1) 

where a and (3 depend on depolarization factor (A) and are detailed in 



Appendix A For instance, one approximation of these parameters are: 



a = l/[A + l/{e\/e e — 1)] and /3 can take two values, we will use: 

£i 

/3sbg = — 3 (2) 
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the other approximation can be found in the Appendix A From now on 



we call this model as M-B. The depolarization factor intrinsically has the 
information of shape and orientation of inclusions. For instance, for elliptical 
inclusions with 90° of orientation respect to the applied electric field, aspect 
ratio (semimajor axis a / semiminor axis b) a/b = 1/3 and E\je e = 20/2, 
Mejdoubi and Brusseau obtained A equal to 0.732 or 0.772, using (3mg or 



Psbg, as defined in Appendix A , respectively. It is remarkable that these 



values were fitted with 0; < 0.1. 

2.3. Micro-structure Parameters 

Trabecular bone structure has been extensively studied and there are 
several parameters to characterize it (see for instance [15]). In this work 
four parameters are studied from 2D images of trabecular bones: BV / TB, 
trabecular thickness, mean intercept length (MIL) and degree of anisotropy 
(DA = MIL max /MIL min ). The former is asociated with bone porosity and 
the latter with anisotropy. For BV / TB and MIL calculus, see [12] and 
references therein. For the trabecular thickness calculus we trace n parallel 
lines in ortogonal direction with respect to the main trabecular direction (we 
define it as angle ?7(MIL max )), measuring the thickness of the intersections 
with trabeculae. 



2.4- Images 

Images used here are avaiable in [29] and in reference [17]. These images 
are used to build the realistic dielectric models and image processing is ap- 
plied in order to define the two-phase material [30 J. 

Samples used here are named as sample 1, 2, 3, 4, 5 and 6, they are all 
shown in section [3] (Figj6|. Figure [2] shows just some of them. Figure [2] (a) 
shows trabecular bone tissue from an old person (sample 1), its porosity is 
approximately 0.75. Polar plot in Fig|2] (d) shows the corresponding MIL 
parameter (r] is the slope angle of the parallel lines for MIL calculus). It 
seems to be isotropic, the maximum of the MIL parameter is around 0° but 
it almost a circle. Figure [2] (b) shows a highly porous old bone as well with 
porosity equal to 0.79 (sample 2). Clearly its MIL maximumm is around 90° 
(see Fig|2] (e)). Figure [2] (c) shows a young bone (sample 3) with porosity 
equal to 0.46 and 0° main direction (Fig|2](f)). 
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Figure 2: Microscopic characterization of main direction: mean intercept length calculus, 
(a), (b) and (c) correspond to sample 1, 2 and 3 as defined in the text (see subsection 2.4), 
respectively, (d), (e) and (f) are the polar plots of the MIL parameter and correspond to 
sample 1, 2 and 3, respectively, rj is the slope angle of the parallel lines for MIL calculus. 
The lines are guides for the eyes. 



3. Results 

3. 1 . Validation 

Validation of method is performed using an hexagonal array of discs, it 
has an analytic solution [31] (see Figj3|. Figures [3] (a) and (b) correspond 
to inclusions of air embedded in a matrix with e e = 16 and a e = 0.043 S/m. 
Figures [3] (c) and (d) correspond to a potential model of trabecular bone tis- 
sue. Inclusions are the trabeculae and the matrix is the bone marrow. These 
validation media were also used by Bonifasi-Lista and Charkaev [4]. Values 
are e e = 5.485 (bone marrow), E\ = 20.584 (trabeculae) and o~ c = 0.043 S/m 
(bone marrow), = 0.364 S/m (trabeculae). The fraction occuped is varied 
by changing the disc radius. A good agreement is observed. 
Disc shape inclusions randomly located are also tested (see Figj4| (a)) with 
similar values to those used in reference [21]. An excelent agreement is ob- 
served. Different circle radii are tested, it can be seen that permittivity and 
conductivity do not depend on the radius but on the fraction occuped ((/>i). 
It is known that not only fraction occuped by the inclusions determines the 
effective dielectric properties but also their shape [22]. Ellipse shape inclu- 
sion media are tested with effective formulaes obtained by reference [22J (see 



Appendix A). The only difference to Mejdoubi and Brosseau work is that, 



here, orientation of the ellipses is a random angle with normal distribution 
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Figure 3: Effective permittivity and conductivity of hexagonal array discs as a function 
of fraction occuped (</>i). Comparison between analytic solution (continuous line) and 
simulation (v). (a) and (b) correspond to the effective permittivity and conductivity 
of a medium with e c = 1.0, 16.0 and ct;, ct = 0.0, 0.043. (c) and (d) correspond to 
the effective permittivity and conductivity of a medium with e\, e c = 20.584, 5.485 and 
o-i, cr c = 0.364, 0.043 



7 







(d).- 







Figure 4: Effective relative permittivity of media with disc and ellipse shape inclusions 
randomly located, (a) Medium with discs for different radius (ri, r 2 and r$, V, □ and 
o, respectively) compared with analytic solution of hexagonal array (continuous line). 
r\ < r 2 < r 3 . (b) Effective permittivity of discs (o) compared to ellipses at 7 = 45° (A) 
and Wiener (continuos and dash-dotted lines) and H-S (cross and dotted lines) bounds. 
(c)-(d) Effective permittivity of oriented ellipses at 7 =0° (V), 90° (□) compared to M-B 
bounds (dotted and continuos marks for /3mg, cross and dashed-dotted marks for /?sbg)- 
In (a), (b) and (c) £i,£ c = 1, 16. In (d) £i,£ c = 16, 1. 



(7 ~ -AT (7, <y 2 = 100)). Three mean values of orientation are evaluated 7 = 
0°, 45° and 90°. For all mentioned above cases overlaping of inclusions is al- 
lowed. Figure [1] (b) shows comparison between ellipses (7 = 45°) with discs 
and effective bounds (Eqs.A.l A.4). In Figures [4] (c) and (d) comparison 



between ellipse shape inclusion media and bounds (EqjTJ is shown. Even 
these models were fitted for 0; < 0.1, a relative good agreement is observed 
for fraction occuped < 0.2 if /?sbg formula is used (see dashed-dot and cross 
marked in FigsQ(c) and (d)). 

Figure [5] (a) shows an example of a model totally anisotropic. Ellipses are 
oriented and arranged in parallel lines. The slope of the lines is varied an 
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(a) (b) (c) 

Figure 5: Ellipse models. Percolated (o) and not percolated (□) models. The lines are 
guides for the eyes, (a) Simulation slice examples, (b)-(c) Polar plots for effective relative 
permittivity and effective conductivity, respectively. 9 is the slope angle of the oriented 
ellipses which is coincident with the slope of the lines (see text for explanation). 



angle 6 (see Figj5] (a)) and the effective dielectric properties are calculated. 
Source is always y polarized, named E y , and ~ 0.054. Two kind of mod- 
els are tested, allowing (percolated) and not (not percolated) overlaping of 
ellipses. Note that, at the simulation frequency of this work (950MHz), per- 
colated and not percolated systems have a similar behaviour. M-B model 
using /?sbg gives depolarization factor A = 0.901 for E y _L MiX max and A = 
0.120 for E y n MIL max , which gives e eS = 5.729, 6.120 and a cS = 0.046, 
0.054 S/m, respectively. Values obtained by numerical simulations are e c s = 
5.689, 6.120 and a e s = 0.045, 0.053 S/m. We see that the agreement between 
simulations and effective models is relatively good. 

3.2. Simulation of realistic models 

Figure [6] shows the samples used here to build the realistic models. Models 
are simulated as a two phase medium where trabeculae are the inclusions and 
the background is the bone marrow. The images are rotated an angle 9 (at 30 
degree steps) and the simulation is performed (the same method, as explained 



in subsection 2.1, is used to calculate the effective dielectric properties). The 
applied field is always y polarized. 

Table [j] resumes their micro-structure parameters (calculated as described in 



subsection 2.3). There is only a young sample (sample 3). Sample 1 and 3 
have the lower DA values. 

Table [2] shows the dielectric properties of samples we have obtained. max y 
#min are the angles for which dielectric properties are maximum (£ e ff max and 
cr e ff max ) and minimum (e e ff min and cr e s min ), respectively. It is remarkable that 
e c s and (j e ff are maxima when the applied field is totally aligned with the main 
direction of the sample (r/(MIL max )). In this table the ratios £ e fT max /£cff min and 
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sample 4 sample 5 sample 6 



Figure 6: Samples. 



Table 1: Microstructure properties. 



Sample 


Porosity 


BV/TV 


^(MIL max ) 


Trab. Thick. [mm]±std 


DA 


1 


0.75 


0.25 


0° 


0.11±0.16 


1.08 


2 


0.79 


0.21 


90° 


0.18±0.13 


2.23 


3 


0.46 


0.54 


0° 


0.23±0.32 


1.14 


4 


0.88 


0.12 


130° 


0.13±0.08 


1.53 


5 


0.74 


0.26 


90° 


0.21±0.19 


1.40 


6 


0.74 


0.26 


90° 


0.13±0.14 


1.42 
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Table 2: Effective dielectric properties obtained by simulation of realistic models. 



Sample 


a 

"max 


^min 


£ eff max 


£pff ■ 
Lii mm 


^Cffmax 
cii mm 




(Tpff . 

Cii mm 


CTpff . 

Lli mm 


1 


90° 


30° 


7.503 


7.096 


1.06 


0.077 


0.0706 


1.09 


2 


0° 


90° 


6.950 


6.247 


1.12 


0.069 


0.054 


1.28 


3 


90° 


0° 


11.446 


10.345 


1.11 


0.154 


0.149 


1.18 


4 


120° 


30° 


6.119 


5.864 


1.04 


0.053 


0.048 


1.10 


5 


0° 


60° 


7.312 


6.798 


1.07 


0.075 


0.063 


1.19 


6 


0" 


90° 


7.189 


6.782 


1.06 


0.071 


0.063 


1.13 




(e) (f) 



Figure 7: Effective dielectric properties of sample 2 (a-c) and sample 4 (d-f ) as a function 
of angle (9). (a) and (d) show the samples and main direction, (b) and (e) effective relative 
permittivity, (c) and (f) effective conductivity. The lines are guides for the eyes. 

cr e ff max /<T e ff min are also shown. It is observed that the latter always gives larger 
values than the former. See for instance: sample 2, it has the maximum 
degree of anisotropy (2.23) and the maximum of conductivity ratio (1-28) as 
well, while permittivity ratio is 1.12. Figure [7] shows tipical simulations of 
sample 2 (Fig|7] (a)) and 4 (Fig|7| (b)). Polar plots shows effective relative 
permittivity and conductivity as a function of 6. (b) and (c) correspond to 
sample 2 and (e) and (f) correspond to sample 4. It should be noted that 
main direction can be detected by both permittivity and conductivity. 

3.3. Effective models 

In this subsection comparison between effective models and realistic mod- 
els are performed. It was shown in [22] that depolarization factor (a scalar 
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Table 3: Effective M-B models using /?sbg- 



a/b 



0.12 



0.21 



BV/TV 
0.25 



0.26 



0.54 



1.0 



1.5 



Ocff 

£ cft'max 

£efF ■ 
mm 

CT cffmax 

^eff ■ 
mm 

£ cff m ax 

£eff ■ 
mm 

^cffmax 



6.688 
0.065 



6.471 
6.245 
0.057 
0.053 

6.744 
6.129 
0.064 
0.051 



1.04 
1.09 

1.10 
1.25 



7.103 
0.066 



7.342 
6.910 
0.071 
0.062 

7.844 
6.689 
0.085 
0.058 



1.06 
1.14 

1.17 
1.47 



7.342 
0.071 



7.565 
0.073 



7.873 
7.315 
0.080 
0.068 

8.512 
7.028 
0.098 
0.062 



1.07 
1.19 

1.21 
1.58 



10.703 
0.124 



in 2D) for media with elliptical inclusions depends on aspect ratio (a/6) and 
orientation. As it is shown in Figj4] effective M-B models have relative good 
agreement for ; < 0.2 when /3sbg is used. In general this result was ob- 
served for several simulations not shown here. Table [3] shows some results 
of M-B models (using /3sbg) as a function of 4>\. We replace 0j by BV/TV 
parameter of each sample. Three aspect ratios are observed: ? =1.0,1.5 and 
3.0. Note that, for high values of DA, although the absolute values agree 
relatively, the correlation with the relative values of | =1.5 is even better 
(values highlighted in boldface). 

4. Discussion 

Results in this paper are limited to a central frequency (950MHz). This 
choice was based on the fact that most of the works on microwave tomography 
use this frequency. Moreover, practical implementation are generally devel- 
oped in frequency domain, based on network analyzer plus antenna array. 
Therefore we assume that constant values of permittivity and conductivity 
of each component of bone tissue (trabecular and bone marrow) are a good 
representation of a real measurement. 

Recent works have begun to study microwave imaging of bone tissue in vitro 
[T0"l[l"2"] and in vivo [11J. In general, the current interest is to find correlations 
between dielectric properties and bone mineral density or bone volume. It 
is surprising that the relationship between anisotropy and dielectric proper- 
ties has been relatively little studied. We will mainly focus this discussion 
into two aspects: firstly, we will present the correlation between BV/TV and 
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effective dielectric properties and, secondly, we will comment the results on 
anisotropy. 

Comparison between dielectric properties and bone volume of Table [T] and 
[2] shows a positive relationship. This means, the higher BV/TV the higher 
effective permittivity and conductivity. This is very intuitive and is in ac- 
cordance with two phase models: when the fraction of inclusion (with higher 
relative permittivity) grows so do the permittivity and conductivity. Differ- 
ences are found when this simulation results are compared to measurements 
in vivo. In reference [11] the authors found that an affected calcaneus (with 
a lower bone volume) has a 23% higher relative permittivity as compared to 
a normal one. It is necessary to remark that it is just one sample without 
statics support. The same group has recently published a paper on in vitro 
measurements of trabecular porcine femur samples (six bone specimens in 
saline solution) [12]. They studied variation of dielectric properties as a func- 
tion of mineralization. Degree of mineralization was altered (not naturally) 
by acid treatment, their results are in agreement with their previous work 
It also agrees with [13], in which degree of mineralization was measured 
and negative correlation between bone mineral density and relative permit- 
tivity was found (six samples of human osteoporotic trabecular samples in 
saline solution, it is also a small number of samples). 

Some simple calculation can be performed with the help of effective models. 
For instance, using Eq|T]for a total aligned electric field and bone orientation 
with 0; = 0.26, £i = 20.584, a { = 0.364 and e e = 78.00, a e = 1.39 (simi- 
lar values to 0.9% saline solution) gives 61.76 and 1.08, for permittivity and 
conductivity, respectively. For fa =0.54 gives 42.79 and 0.749. These results 
give a possible explanation of the relation between BV/TS and dielectric 
effective properties (at least for in vitro measurements): samples with lower 
BV/TV have more space for water or saline solution, which increases their 
permittivity and conductivity. Special care must be taken for interpreting in 
vivo measurements. What is taken the place of bone in a resorption process: 
Is it predominantly some high water content connective tissue/cells/blood 
(high dielectric properties) or is it bone marrow (low dielectric properties)? 
As first image in vivo has shown [TT]. it seems that a loss of bone tissue in- 
creases both permittivity and conductivity, then first option could be the 
more feasible process. 

Regarding anisotropy, Table [2] shows that polarized electric field totally 
aligned with main direction of the sample gives the higher relative permit- 
tivity and conductivity. On the contrary, minimum values are found if the 
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electric field and main direction are ortogonal. 

Table [2] also shows that degree of anisotropy correlates relatively well with 
permittivity ratio (s c s m ^/£ c s min )- Maximum value of DA (2.23) for sample 2 
gives also maximum ratio of relative permittivity (1.12). In this case, per- 
mittivity is 12% higher when the electric field is totally aligned with the 
main direction. Conductivity ratio (rigth most column of Table [2J shows the 
strongest relationship with DA, in the case mentioned before conductivity is 
28% higher when there is a total alignment. As a conluding remark of these 
results: conductivity could be a better estimate of DA. 
Effective ellipses models are more approximated when fraction occuped is 
low, strictly speaking for 0; < 0.1. Table [3] resumes the application to tra- 
becular bone tissue. It can be seen a relative good agreement for BV/TV < 
0.26 (see Table [3] values in boldface). For instance, sample 2, 4, 5 and 6 have 
similar values to those given for a ellipses model with aspect ratio a/b = 
1.5. It shows an interesting trend: ellipse models represent well samples 
with high degree of anisotropy and low fraction occuped. On the other hand, 
when samples have low DA and low fraction occuped circular inclusions mod- 
els seem to be better. 

As a final conclusion, we can say that samples with high degree of anisotropy 
can give until 20% of difference between aligned and ortogonal conductivity 
values. Relative permittivity shows lower variation. Semenov et al. [TU] 
studied microwave-tomographic imaging of the high dielectric-contrast ob- 
jects (as bone tissue and muscle or skin). Several difficulties of reconstruc- 
tion algorithms were found by the authors even by taking bone as isotropic 
tissue. Microwave imaging of bone tissue is still in its begining, but according 
to the results presented here, anisotropy must be taken into account when 
measuring trabecular bone tissue. 
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Appendix A. Effective approaches 

Wiener approach (W) 

£efr max = + (1 - <Pi)e e (A.l) 
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£pff 



[1 - 0i)£i + 



Hashin-Strikman approach (H-S) 



£i-e 2e e 



£ P ff . = £\ H 7 

eii mm 1 ] 



l-0i 



+ P~ 



Mejdoubi-Brosseau approach (M-B) 



£ e ff = • / ( — , <t>h A ) « £ e (1 + «0i + 



(A.2) 

(A.3) 
(A.4) 

(A.5) 



with a = 1/L4 + 1/ (si/e c — 1)] and /3 can take two values depending on the 
approach: 



or 



P. 



MG 



2 1 

A + T. + -7 



3SBG 



A 2 I — - 1 



A[ --1 



V 

/ 
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